Domain wall propagation through spin wave emission 



(N 

o 

C 



03 



o 
o 



> 

O 

(N 



X 



X.S. Wang 1 , P. Yan 2 , Y.H. Shen 1 , G.E.W. Bauer 3 ' 2 and X.R. Wang 4 >£] 

1 Physics Department, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong 
2 Kavli Institute of NanoScience, Delft University of Technology, Delft, The Netherlands 
3 Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan and 
4 School of Physics, Wuhan University, Wuhan, P. R. China 

We theoretically study field-induced domain wall (DW) motion in an electrically insulating fer- 
romagnet with hard- and easy-axis anisotropies. DWs can propagate along a dissipationless wire 
through spin wave emission locked into the known soliton velocity at low fields. In the presence 
of damping, the mode appears before the Walker breakdown field for strong out-of-plane magnetic 
anisotropy, and the usual Walker rigid-body propagation mode becomes unstable when the field is 
between the maximal-DW-speed field and Walker breakdown field. 
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Magnetic domain- wall (DW) propagation in nanowires 
has attracted attention because of the academic interest 
of a unique non-linear system [H-Ql and potential ap- 
plications in data storage and logic devices [JQ- The 
field-driven DW dynamics is governed by the Landau- 
Lifshitz-Gilbert (LLG) equation which has analytical 
solutions in limiting cases [l|, @] , such as the soliton solu- 
tion Q in the absence of both dissipations and external 
magnetic fields. The interplay between spin waves (SWs) 
and DWs has also received attention, including DWprop- 
agation driven by externally generated SWs [9l4lll| and 




SW generation by a moving DW [12|, |13(. Our under- 
standing of the field-induced DW motion is nevertheless 
far from complete. According to conventional wisdom 
DWs move under a static magnetic field only in the pres- 
ence of energy dissipation [J Numerical evidence 
against this view therefore came as a surprise [l2j ]. 

We report here a physical picture for the SW emission- 
induced domain wall motion for a head-to-head DW in 
a magnetic nanowire with easy axis along the wire (z- 
direction) as shown in Fig. Q] Let Kn and K± be 
anisotropy coefficients of the easy and hard axis (along 
the rr-direction) , respectively An external field along 
the wire rotates the DW out of the j/z-plane. The DW 
structure thereby experiences an internal field in the x- 
direction twisting the DW plane and generating a non- 
uniform internal field along the wire. This field causes pe- 
riodic deformations of the DW structure, such as "breath- 
ing" [l| by which the entire DW precesses around the 
wire axis while its width shrinks and expands periodi- 
cally. The local modulation of the magnetization texture 
generates SWs (wavy lines with arrows in Fig. 1) that 
radiate away from the DW center. The energy needed to 
generate the SWs has to come from the Zeeman energy 
14| that is released by propagating the DW. The DW ve- 
locity of a dissipationless ferromagnct in the steady state 
may then be expected to be proportional to the SW emis- 
sion rate. 

In this Letter, we numerically solve the LLG equation, 
initially without damping in order to confirm the above 



FIG. 1: (color online) Schematic of a head-to- head DW of 
width A in a magnetic nanowire. H is an external field along 
wire-axis defined as ^-direction. DW breathing and other 
types of periodic DW texture deformations emits spin waves, 
denoted by the wavy lines with arrows. 



mentioned relation between spin wave emission and DW 
propagation. Depending on K± and the magnetic field, 
breathing or more complicated periodic transformations 
of the DW emit spin waves. The DW propagation at 
low fields tends to lock into a particular soliton mode in 
which the energy dissipation rate due to the SW emis- 
sion is balanced by the Zeeman energy gain. We predict 
robust spin wave emission that persists in the presence of 
Gilbert damping and renders the usual Walker rigid-body 
propagation mode unstable in region below the Walker 
breakdown field. 

The LLG equation reads 



dm 

— - = m x heff 

at 



dm 

am x — , (I) 



where m is the unit direction of the local magnetiza- 
tion M = mA/ s with saturation magnetization M s and 
a is the Gilbert damping constant. The effective mag- 
netic field of our biaxial wire (see Fig. I) is h off = 
K\\m z z — K±m x x + Ad 2 m/dz 2 + Hi, consisting of in- 
ternal and external fields in the unit of M s . A is the 
exchange constant. Time, length and energy density are 
measured in units of (7-M s ) with gyromagnetic ratio 
7, the DW width at equilibrium Ao = nWA/Kii , and 
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FIG. 2: (color online) Snapshots of m± = ^/m| + for 
= 0.01 at various times t for (a) K± = 4 and a — 0; (b) 
A' ± = 10 and a = 0; (c) K± = 4 and a = 0.001; (d) A ± = 10 
and a = 0.001. 



FIG. 3: (color online) Ax-dependence of the DW velocity v 
(in units of tvjM s ^/A/K^) for fixed H = 0.01 and a = 0. 
Insets: Time dependence of the DW center position (dashed 
line) and the azimuthal angle <j> (solid line) of magnetization 
at the DW center for K± = 4 (left) and K± = 10 (right). 



/i Ai^, respectively. We chose parameters of the elec 



trie insulator Yttrium Iron Garnet (YIG) with [111, |15| 
A = 3.84 x 1(T 12 J/m, K\\ = 2 x 10 3 J/m 3 , 7 = 35.1 
kHz / (A/m) and M s = 1.94 x 10 5 A/m, and the cor- 
responding time and length units are 1.46 x 1CP 10 s and 
1.38 x 10 -7 m, respectively. K± and a are treated as ad- 
justable parameters depending on the sample shape and 
microscopic order. We solve Eq. ([T} by a numerically 
stable method 16j. The mesh size is chosen to be 0.009, 



corresponding to the YIG lattice constant (1.24 nm) 

To prove that a DW under an external field indeed 
emits spin waves, we plot the snapshots of the distribu- 
tion of mj_ = ^/to 2 + Trig for K± = 4 and 10 (in units 

of /i(iM s 2 ), a = 0, and H = 0.01 (on at t = 0) at various 
times in Figs. 2a and 2b. At t = 0, right before the 
external field is switched on, m± follows the Walker DW 
profile [l[, the DW center is located at the center (z = 0) 
of a wire with length I8OA0, while the DW magnetization 
lies in j/z-plane. Curves are offset for better visibility. As 
time proceeds (t > 0), SWs (wavy features) are emitted 
into both directions, while the DW center (peak) moves 
simultaneously along the field slower than the SWs. The 
velocities v for a fixed magnetic field increase monoton- 
ically with K± (Fig. [3]). It does not depend sensitively 
on small K± (< 4), but grows rapidly when K± is close 
to 8, and becomes an almost linear function of field for 
K± > 12. From the time-dependence of the position 
(dashed lines) of DW center and its azimuthal angle (solid 
lines) shown in the insets of Fig. [3] for a small and large 
K±, we trace the periodic DW deformations in the dif- 
ferent field regions. We note that such large K± value 
may be realized in YIG samples subject to mechanical 
strains [17|. 

At a small K± = 4 (left inset), m at the DW center ro- 



tates around the wire, while the center position oscillates 
back and forth but also moves slowly along the applied 
field. This oscillatory motion is synchronized with the 
breathing, i.e., the periodic energy absorbing and releas- 
ing in the form of periodic oscillation of the DW width 
A = ■K^/Ajjjq+Yl cos 2 4>) fl, where </> is the tilt an- 
gle of the DW plane (with equilibrium value tt/2). This 
breathing excites spin waves as shown in Fig. 2a. Low 
velocities corresponds to weak SW emission. At a large 
K± = 10 (right inset), the azimuthal angle of the DW 
center approaches a fixed value and the DW center po- 
sition moves at a constant velocity since 4> and the DW 
energy are almost constant. In this case, m still rotates 
around the wire axis while the DW center propagates 
along the wire with a fixed <fi. The large K± twists the 
DW plane into a chiral screw-like structure that changes 
shape periodically during the magnetization precession 
while the DW center "drills" forward. This drilling mode 
is much more efficient in emitting SWs than the breath- 
ing mode, leading to a relatively high DW propagating 
speed. For YIG parameters the SW velocity exceeds that 
of the DW, therefore, in contrast to Ref. 12( , we observe 
bow as well as stern SW excitations. 

Figure [4^, displays the steady state DW velocity in a 
dissipationless wire with K± = 4, 10, and 16 with pa- 
rameters otherwise identical to those of Fig. 3. As a 
function of field it increases abruptly for small values, 
reaches a maximal value, and decreases again. When 
K± is reduced from 16 to 4 the DW changes from a 
drilling to the breathing motion, resulting in a signifi- 
cant drop of DW velocity. The decrease of DW veloc- 
ity with field should not be interpreted as suppression 
of SW emission or damping of DW propagation by spin 
wave emission (l2l [l8| . The Zeeman energy released by 
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the DW motion at a rate 2HMv [14| should be equal 
to the energy rate carried away by the SWs. Therefore, 
provided that the latter increases sub-linearly with H 
(oc H 13 and /3 < 1), the DW propagation speed must 
decrease with field. The initial rapid rise of the DW 
velocity at small fields is related to the soliton solution 
lntan(0/2) = tt(z - vt)/A of the LLG Eq. (1) with soli- 
ton velocity v = —\JA/ (K\\ + K± cos' 2 4>)K± sin 2<j) for 
a = and H = [8[, where 9 is the polar angle of m. 
This can be seen from the plot of — sin 20 for the sat- 
urated (f> vs. v-s/jKn + K± cos 2 <fi)/A/K± (symbols) for 
H G [0, 0.0009] in the inset of Fig. Ha for K± = 16. 
The numerical data agree precisely with the soliton ve- 
locity formula (solid line) . Since conventional solitons do 
not satisfy the LLG equation in the presence of an exter- 
nal field we unearthed here a hybridized mode of solitons 
and spin waves. In contrast to the zero field case, there is 
only one particular soliton mode in which the SW emis- 
sion is balanced by the Zeeman energy change, viz. the 
moving DW. This holds as long as the field is smaller 
than the value at which — sin 2<j) reaches its maximum 
value 1. Beyond that field, the soliton mode becomes in- 
stable since the SW emission rate cannot keep up with 
the released Zeeman energy and other propagation modes 
have to take over. This soliton instability point is em- 
phasized by vertical bars in the data points for K± = 4 
and K± = 10, indicating the threshold fields below which 
the soliton formula holds. For K± = 16 this value is out 
of the range of Fig. [4^. 

For YIG parameters and K± = 10, the DW velocity 
at H = 0.01 (~ 24 Oe), is about 140 m / s. For compari- 
son, the corresponding DW velocity by Walker rigid body 
propagation for the same parameters and a = 0.05 is 57 
m / s [l| . We may conclude that the DW velocity in low 
dissipation magnetic insulators is of the same order of 
magnitude as in highly dissipativc ferromagnetic metals. 

DW propagation through spin wave emission exists in 
any magnetic wire with transverse magnetic anisotropy 
irrespective of the Gilbert damping. Figures ^ and [5Ji 
look very similar to Figs. [2jt and[2jD in spite of the finite 
damping a = 0.001. Naturally, in the presence of damp- 
ing the SWs can propagate only over finite distances, 
which explains why they have been overlooked in most 
experimental and numerical studies. In Fig. 2b (or 2a) 
and Fig. 2d (or 2c), the DW velocity is higher in the 
presence of a small non-zero damping, since damping dis- 
sipates an energy on the top of the SW emission and the 
DW velocity is proportional to the energy dissipation rate 
141 ] . Here we find a mixed DW propagation mode that 
profits from both Gilbert damping and spin wave emis- 
sion. Figure|4jD is the field dependence of the DW velocity 
for a = 0.001 and various Kj_, while other model param- 
eters are unmodified from Fig. |4ji. Except at very small 
fields, Fig. 2}} is similar with Fig. The numerically 
obtained velocities (symbols) in Fig. 0J) agree well with 
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FIG. 4: (color online) (a)(b) Field dependence of DW veloc- 
ity for various K±. Field and velocity are in units of 0.01Af s 
and n^/ A/ K^jM s , respectively. A, K» and M B are YIG pa- 
rameters, (a) a = 0. Inset: symbols display — sin2</> vs. 
vy/(K\\ + K± cos 2 <f>) I A /K± , where v is the soliton velocity. 
The line is the linear relationship from soliton theory, (b) 
a = 0.001. The solid curves are the corresponding Walker 
solutions, (c) Snapshot of mi for H — 0.007 and K± — 16 
at t — 400 (for the datapoint indicated by a circle in (b)). 
The numerical results deviates from the Walker profile that 
is clearly narrower. Inset: Expanded view of m± close to the 
DW center. The relative deviation from the Walker profile 
(solid curve) is pronounced and spin wave emission is con- 
spicuous. 



Walker's rigid-body propagation (solid curves) [Jill] be- 
low some maximum field Hm — 2(K\\/K±)°- 25 Hw for 
Kx » K\\ i where Hw = aK±/2 [l[ is the Walker break- 
down field. However, deviations are obvious for fields 
between Hm and Hw (where the solid lines end). This 
implies that Walker rigid-body propagation is not stable 
for H G [Hm, Hw] with respect to the spin wave emission 
mode. 

In order to prove numerically that the Walker solution 
is only stable for H < Hm, we solve the LLG Eq. ([1]) 
starting from an initial magnetization configuration that 
deviates slightly from the rigid-body propagation mode. 
The magnetization indeed converges to the Walker profile 
for H < Hm- However, this changes when we consider 
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H = 0.007 for K± = 16 (indicated by dashed circle in 
Fig. |4)d) with all other parameters the same. According 
to Walker theory [lj], H M = 0.00344 and H w = 0.008 
(indicated by dashed lines in Fig. The symbols in 

Fig. |3t give a snapshot of m± for t = 400, at which 
the transients die out and the DW center propagates to 
about z = 97. The distribution deviates significantly 
from the Walker profile (solid curve). The spin wave 
emission is clearly observed in the expanded tail in the 
inset of Fig. Hfc. Exact analytic solutions do not make 
numerical methods obsolete. 

There are several corollaries of the DW propagation 
mode by spin wave emission. The emitted spin waves 
from one DW, for example, can mediate an attractive 
force on the nearby DW, since a DW moves against pass- 



ing spin waves by spin transfer jllj . This causes crosstalk 



in wires with more than one domain walls with conse- 
quence for the "race track" memory 0]. This DW-DW 
attractive force has a finite range governed by the Gilbert 
damping and is sensitive to material parameters and ge- 
ometry. The increase of the effective damping by SW 
emission [l~8[ is not restricted to DWs but will appear 
in any time-dependent magnetization texture including 
magnetic vortices and cannot be captured by the Gilbert 
phenomenology. On the other hand, our results possibly 
open alternatives to manipulate and control the effective 
damping in magnetic nanostructures. 

In conclusion, we prove that a DW in a wire with a 
finite transverse magnetic anisotropy undergoes a peri- 
odic transformation under an external magnetic field that 
excites SWs. The energy carried away must be compen- 
sated by the Zeeman energy that is released by DW prop- 
agation along the wire. The DW propagation adopts one 
particular soliton velocity at low fields. This SW assisted 
DW propagation can be attributed to the SW emission 
generated by DW breathing and drilling modes. In the 
presence of damping it competes with and appears before 
the Walker breakdown. Also, the spin wave emission will 
mediate a dynamic attractive force between DWs with a 
range controlled by the Gilbert damping. 
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